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Abstract 

We present a computer simulation of exciton-exciton scattering in a quantum well. Specifically, we use quantum Monte 
Carlo techniques to study the bound and continuum states of two excitons in a 10 nm wide GaAs/Alo.3Gao.7As quantum 
well. From these bound and continuum states we extract the momentum-dependent phase shifts for s-wave scattering. A 
surprising finding of this work is that a commonly studied effective-mass model for excitons in a 10 nm quantum well 
actually supports two bound biexciton states. The second, weakly bound state may dramatically enhance exciton-exciton 
interactions. We also fit our results to a hard-disk model and indicate directions for future work. 
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1. Introduction 

Exciton-exciton interactions in quantum wells are 
becoming increasingly important to science and tech- 
nology. For example, the interaction of two excitons 
is a mechanism for non-linear optical response [1]. It 
is the excitonic component of polaritons that lead to 
their interactions. Also, in experimental studies of cold 
excitonic gases in quantum wells, exciton-exciton scat- 
tering is a process that drives the excitons towards a 
Bose distribution. 

Exciton-exciton interactions are difficult to theoret- 
ically predict [2] . Fast collisions between excitons may 
be treated with time dependent perturbation theory 
[3] . Carrier exchange is known to be a significant con- 
tribution to the scattering process [3]. At a different 
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extreme, slow moving excitons interact through dipole 
interactions. If the excitons have intrinsic dipoles, 
perhaps from carefully engineered quantum well band- 
structure or an applied electric field, the repulsive 
dipole-dipole force can dominate interactions [1,4]. In 
the more common case of unpolarized excitons, quan- 
tum fluctuations of the polarization leads to attractive 
van der Waals interactions, which are often neglected. 
Finally, we know that the wavefunction describing 
elastic exciton collisions must be orthogonal to any 
bound biexciton states. The presence of weakly bound 
biexciton states can have dramatic effect on exciton 
scattering cross-sections [2]. 

In this work we have extended a quantum Monte 
Carlo technique in order to to study exciton scatter- 
ing in a GaAs/Alo.3Gao.7As quantum wells. An ear- 
lier paper studied bulk exciton-exciton scattering for 
a generic model of a semiconductor [2]. The simula- 
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tion technique directly samples the energy of the four- 
particle scattering wavefunction. This allows all mech- 
anisms listed in the previous paragraph — carrier ex- 
change, dipole and van der Waals forces, and orthog- 
onality to bound biexcitons — to be fully incorpo- 
rated in the numerical simulation. Random walks are 
used to project essentially exact energies from a set 
of trial wavefunctions, thus providing numerical esti- 
mates of the phase shifts and scattering cross sections 
for elastic exciton-exciton collisions. Our new calcula- 
tions apply this technique to a quantum well, allowing 
a detailed study of the exciton-exciton collision pro- 
cess within the effective mass approximation. Two key 
results of these calculations are the observation of a 
second, weakly-bound biexciton and the numerical es- 
timate of a large cross-section a > 100 nm for low- 
energy exciton-exciton scattering in a quantum well. 

In the next section of this paper, we describe the ef- 
fective mass model we use. After that we give a brief 
description of our computational technique and details 
of our simulation of a GaAs/Alo.3Gao.7As well. The 
last section is a discussion that emphasizes new results 
of this paper, some subtleties of two-dimensional scat- 
tering, and future directions for this research. 



2. Model 

We use a single-band effective mass model, 
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Here the electrons are labeled 1 and 2 and the holes 
are labeled a and b. The isotropic effective masses 
are ml — 0.0667mo and m£ = 0.34mo, and the di- 
electric constant is e = 12. (Future calculations will 
include the anisotropy of the heavy-hole mass.) The 
confining potentials V e and 14 are finite square wells 
with offsets AV e = 216 eV and A 14 = 133 meV. 
These material constants have been chosen to match 
other theoretical studies of excitons and biexcitons in 
GaAs/Alo.3Gao.7As quantum wells for direct compar- 
ison with this work [5,6]. 
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Fig. 1. Projection of eigenvalues with a random walk in 
imaginary time (shown for the simulation with r cu t = 180 
nm). The projection over an imaginary time 0.06 mev - 1 ■ h 
takes a buffer of 200 Monte Carlo steps. 

The electron effective mass and dielectric constant 
define natural length and energy scales for the prob- 
lem: the donor effective Bohr radius, a* B = h 2 e/mle 2 = 
9.97 nm, and the donor effective Rydberg, Ryd* = 
h/2m, 1 a* B = 5.79 meV. In this paper we present nu- 
merical answers to three digits for comparison to other 
theoretical papers; such accuracy is beyond the limits 
of the effective mass model and the leading digit should 
suffice for comparison to experiments. 

For a well of width 9.97 nm, we find an exciton bind- 
ing energy of 9.22 meV, using QMC methods described 
in the next section. Previously published QMC calcu- 
lations report energies of about 10.0 meV for the ex- 
citon binding energy at this well width [5,6]. We have 
also calculated the biexciton binding energy, which is 
around 1.87±0.01 meV, in very good agreement with 
experimental [7] and theoretical [8] values. 

We also see a second, very weakly-bound excited 
biexciton state with binding energy 0.032±0.007 meV. 
This state is not likely to be directly observable, but it 
can have a big influence on the scattering length. Also, 
note that this calculation is within the effective-mass 
approximation with isotropic effective masses, and the 
bound state may be an artifact of this approximation. 
Nevertheless, this results suggests that a second bound 
biexciton state is possible in a GaAs/AlGaAs quantum 
well. 
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3. Calculation Technique and Results 

We use random walks to project a set of wavefunc- 
tions in imaginary time, as described in Ref. [2]. We 
first construct variational wavefunctions correspond- 
ing to a bound state and excited state of scattering. 
The excitons can be paired two ways: electron 1 can be 
paired with hole a and electron 1 with hole b, which we 
denote la, 2b, or we can have the pairing lb, 2a. Eigen- 
states are symmetric or antisymmetric combinations, 



<J>±± = $la 



(2) 



These wavefunctions can be related to spin states 
through the symmetry of the wavefunction. In these 
calculations we are just studying the $++ wavefunc- 
tion, since this is the one that has the bound biexciton. 
This is one channel for scattering, and to get cross 
sections for different exciton spin states we will need 
to study also. 

Next we discretize the scattering states by restrict- 
ing the excitons to never be separated by a distance 
greater than R„ . We choose values of R n between 160 
nm and 280 nm. This turns the scattering states into 
standing waves we can study with QMC. Using QMC, 
we project the wavefunction with a random walk. The 
random walk is repeated projection of the wavefunc- 
tion with exp —tH. In Fig. 1 we show the projection of 
eigenvalues from the simulation with i? n =200 nm. We 
clearly see the bound biexciton projected to an energy 
of about -1.9 meV. There is also a much more subtle 
feature: a second weakly bound biexciton state. This 
can be thought of as a quantized vibrational excitation 
in the radial exciton-exciton separation coordinate. 

In Fig. 2 we show the energy eigenvalues for a series 
of different R n values. We see that the two bound states 
are relatively insensitive to the boundary condition, 
while the scattering states systematically decrease in 
energy with increasing R n . These energy curves hold 
all the information about low energy s-wave scattering. 

Since we know the energies of the scattering states, 
we know the relative momentum k. Because we know 
the scattering function has a node at R n , we can infer 
the phase shift, 8(k). We plot the phase shift in Fig. 3. 
Again, the weakly bound vibrationally excited biexci- 
ton shows up, this time as the phase tends to 2n, in- 
dicating two bound states. The scatter in this data is 
the Monte Carlo noise from our simulations. 
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Fig. 2. Energy vs. cutoff. Small, random fluctuations are 
the statistical noise in the Monte Carlo simulations. 
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Fig. 3. Phase shifts calculated from data points in Fig. 2. 

From the phase shift, we calculate the s-wave scat- 
tering cross sections, given by a — 4sin 2 ((5)/fc. Fig. 4 
shows the scattering cross sections from our raw Monte 
Carlo data. At low energy, the cross section becomes 
quite large, a > 100 nm. Two effects may contribute: 
(1) the weakly-bound excited biexciton state at ~ 0.03 
meV can enhance scattering, and (2) two dimensional 
systems have a logarithmic divergence of the scatter- 
ing cross-section at low energy. The statistical error in 
these simulations prevents us from distinguishing be- 
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Fig. 4. Two-dimesional scattering cross sections (with di- 
mension of length) from data points in Fig. 3. 

tween these contributions. 



if the electron and hole spins are each in singlets as 
required by our <3?++ wavefunction. 

Now we must consider spin. We have only simulated 
the <3?++ scattering states. In bulk, three-dimensional 
calculations, the <!> have a much smaller contribu- 
tion to the scattering, since there is now bound biex- 
cition when electrons and holes are spin polarized [2]. 
Using Table III of Ref. [2] and neglecting the $ con- 
tribution, we approximate the scattering cross section 
of paraexcitons as jg the symmetric value, or 1.7 nm. 
Two orthoexictons in a relative S = state have a 
scattering cross section that is \ the symmetric value, 
or 13.5 nm. Recall that this hard disk analysis neglects 
resonant effects from the weakly bound excited biexci- 
ton. 

Future work will include a study of $++ and $ 

scattering channels for different geometries. We will 
also attempt to extract two parameters, an effective 
range and a scattering strength, from our calcula- 
tions. Such a parameterization could provide input for 
models of thermal equilibration rates and properties a 
weakly-interacting Bose condensate of excitons. 



4. Discussion 



Scattering in two dimensions differs significantly 
from three-dimensional scattering [9]. In particular, 
the scattering length can no longer be simply ob- 
tained from an effective range expansion [10]. In three 
dimensions, the scattering length characterizes the in- 
teraction strength for low energy collisions. The corre- 
sponding quantity in two dimensions is dimensionless: 
the scattering is characterized by an effective range 
and a dimensionless parameter. As a simple analysis, 
we will extract just one parameter from our results by 
fitting our phase shifts to a hard-disk model. This anal- 
ysis cannot capture resonant effects from the weakly 
bound state, but it does provide a useful length scale. 

The phase shift for s-wave scattering for hard disks 
of radius a in two dimensions is 

| cot 5 = In +7+ ia 2 fc 2 + C9(fc 4 ), (3) 

where 7 = 0.5772 ... is the Euler-Mascheroni constant. 
A crude fit of our calculated values of -| cot 5 to the 
hard disk formula gives d ~ 2.8a* B with a statistical 
uncertainty of about 10%. This suggests the excitons 
can be approximated as hard disks with radii of 27 nm, 
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